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Magnetic resonance imaging method 

i 

The present invention relates to a magnetic resonance imaging method 
produce successive magnetic resonance images. ! \ 



to 



for 



ca 
so 



MR imaging has been established as a viable and unique imaging technique 
cardiological examinations. Its utility is related to the ability to provide both, precise 
morphology and assessment of global and regional ventricular function. For 
functional assessment, of the heart, steady-state free-precession (SSFP) has 
become an indispensable sequence due to its superb blood-to-myocardium 
contrast (see e.g. Oppelt A, et ai. F1SP: a new fast MRl sequence. Electrbmed 
1986; 54: 15-18.). Compared to spoiled steady-state sequences, SSFP a 
offers improved sighal-to~noise ratio (SNRJ, and it allows for a shorter rhifiimum 
repetition time (TR). However, the advantages of SSFP come at the expense of 
is increased susceptibility to experimental non-idealities. As a result, it is only With 
recent advances in scanner hardware and sequence design that SSFP has 
become widely available. One prominent artifact in SSFP imaging is the sojcajlad 
banding artifact. This results from unaccounted dephasing of tine magnetisation, 
such as from inhomogeneity of the main magnetic field. To reduce susceptibility to 
20 banding artifacts, it is critical to keep TR as short as possible, so that there is less 
time for unaccounted. dephasing to occur. This, in turn, limits the use of j echo- 
planar readout and, accordingly, acquisition efficiency of SSFP imaging isj lower 
than in multiplanar image formation (see Mansfield P. Multiplanar image formation 
using NMR spin echoes. J Phys C 1977; 10: L55-58.). 

25 * 

SSFP. can be applied in a real-time mode or in a triggered mode. In real-time 
mode, data are acquired in a continuous fashion. The lower acquisition efficiency 

of SSFP limits the achievable spatial and temporal resolutions. In the triggered 

• I I 

mode, the idea is to "freeze" the effects of anatomical motion, which would deduce 
30 image quality. This is achieved using cardiac triggering - the coupling of image 
acquisition to a trigger signal - for cardiac studies, generally to an ECGj signal. 
Using the assumption that most biological motion is rhythmic, triggering permits 
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dato acqurred at uHferent cardiac cyc.e S out at the same point ..thin tha cycle 
be combined to improve signal-to-noise ratio in the fine, image or to a gl 
sufficient data to reconstruct a fully-sampled image. Since the images 
mufcpte cardiaccycies are acquired at the same point within the cycle, thele,. 
of motion are minimized. However, residual motion does exist, due to cycle 
cycle variation, and It is manifested as Image artifacts. Therefore, it is desirable 
reduce the overall scan, duration in order to minimize such variability. A series 
images acquired at different points of the cardiac cycle can also be played bad 
cnematic loop fashion, called cine imaging, to provide physicians a way of viewi 
dynamic activities, Triggered acquisitions are generally performed in conjunct.,, 
with a fast imaging pulse sequence for improved temporal resolution Thi lower 
acquisition efficiency of SSFP mentioned previously leads to a longer overa'll sin 
duration, compared to echo-planar imaging sequences, for example. ! T 

i 

An important application of SSFP imaging is in obtaining volumetric time-reUv^d 
images of the heart. For example, multi-slice, multi-breathhold imaging based (an 
SSFP is currently the method of choice for obtaining cine views of the heart] Frdm 
these images, functional parameters can be determined with high accuracy aid 
reproducibility (see e.g. W. U etal. MR assessment of left ventricular function: 
quantitative comparison of fast imaging employing steady-state acduisftil-n 
(FIESTA) with fast gradient echo cine technique. J. Magn. Reson. Imaging, 20 0*2; 
16, p. 559-564). There are. however, limitations associated with multi-slice, muti- 
braathhold acquisftibns. inconsistencies in the multiple breath-holds may lead to 
misregistration of slices. Furthermore, the low bandwidth-time product, if the 
excitation pulses used for short TR SSFP imaging results in large variations of 
excitation angles within the slice. This not only affects the signal-to-noisl and 
contrast-to-noise ratios, it also hampers the reconstruction of multiplanar 
reformats through the heart due to the slab boundary artifact 



To overcome the aforementioned shortcomings, SSFP with true 3D imaging 
previously developed (see B.A. Jung et. al. Single-breathhold 3D-true Flip 
cardiac imaging. Magn. Reson. Med. 2002; 48, p. 921-925). Cine 3D SSFF? of the 
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heart was shown to offer superior signal-to-noise and contrastrto-noise ralios 
compared to mulHslice sequences. However, in that work, compromises bkwien 
image quality, spatial and temporal resolutions had to be made to accomrhodke 
sufficient volume coverage into a single breath-hold. Using an elliptical shytterin 
ky-k* space, cine volumetric data sets may be obtained at relatively high ispa jal 
resolution in a single breathhold. Still, temporal sampling does not satiny I he 
requirements for accurate left ventricular volume calculations. Generally, it leist 
10 temporal frequency harmonics are considered sufficient for accurate Volume 
curve calculation, requiring a sampling rate of at least 20xHR/60 Hz, wkh ijiR 
being the actual heart rate. 



It is an object of the present invention to speed up SSFP dynamic imagirlg. 
particularly by. a high acceleration factor. The acceleration factor is defined jas tie 
ratio between the acquisition time of a conventional scan with full sampling! of tpe 
is data to the acquisition time of an accelerated scan with reduced data sampling. 

This object is achieved by an MR imaging method with the features of claiiVi 1. 
practical use of this invention, as described in detail below, is in accelerating 
cardiac cine 3D imaging in order to obtain a sequence of cine images in ajmu ;h 
shorter time or with a higher temporal resolution for a fixed acquisition duration. 
Another promising use of this invention is in accelerating free-brcsgthir g, 
untriggered (referred to as real-time hereafter) imaging to achieve higher &patal 
and/or temporal resolutions. ; 



25 The present invention is based on the knowledge that the high intrinsic sighal-to 
noise and cohtrast-to-noise ratios of SSFP imaging make It a good candidate 
reduced data acquisition schemes. Moreover, sparse sampling of /c-spaci 
allow shortening of scan time. At the same time, shorter scan time have the 
benefit of reducing the overall power deposited in the subject, which is a cdj 

30 for cardiac SSFP imaging especially at higher field strengths. 



for 
muy 
added 



Jncem 



IBTp002 



i- — r .L.no/m/oni-R iqtif; 



Pmrvf rvr 'K7« P 01 1 



02/07/2003 18:01 



+41-5S-4859495 



SPIERENBURG PARTNER 



10 



15 



20 



25 



30 



4 



However, accelerated SSFP is pro ne to image artifacts resulting from Instable* 
of the gradrent hardware system, such as eddy currents. The gist of this Invention 
is to modify the profile order (also known as the phase-encode order) to reduce 
the said artifacts. These artifacts are often present in SSFP, but they generLy 
remain at low level in practice. Therefore, linear profile order sampling in k-splce 
prevails, particularly in 2D imaging. In the case of high acceleration factors! which 
are desirable in terms of reduced scan time or increased spatlotemporal 
resolutions, these artifacts are amplified drastically. In those cases, the resulting 
image artifacts overwhelm the image contents and must be suppressed. 

Further advantages of the invention are disclosed in the dependent claims land 
the following description in which an exemplified embodiment of the invention 
described with respect to the accompanying drawings. It shows 



Fig. 1 
Fig. 2 



Fig. 3 
Fig. 4 
Fig. 5 
Fig. 6 

Fig. 7 

Fig. 8 

Fig. 8 



regular sampling pattern of magnetic resonance signals, 
time series of fully sampled images 1 and corresponding 'signal 
distribution of object intensities 3 along the dashed line 2 after Inverse 
Fourier transform overtime, 
schematics of the k-t BLAST method, 

Cine 3D SSFP: influence of sampling order on image artifacts, ! 
ReaJ-time 2D SSFP: influence of sampling order on image artifacts 
the sampling pattern for cine 3D k-t BLAST with integrated 
data, 

a cine 3D single breath-hold data set obtained with 5x k-t BLAST in s 
volunteer, 

a cine 3D. single breath-hold data set obtained with 5x k-t BLAST in 
cardiac patient and 

diagrammatically a magnetic resonance imaging system. 



The MR method according to the present invention is based on an approach 
accelerated data acquisition in dynamic imaging, called k-t BLAST (Broad-i 
Linear Acquisition Speed-up Technique) and k-t SENSE (SENSftivity Encoding) 
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The /c-f BLAST and k-t SENSE methods are described in an abstract by -J. Tsao 
et al. (Prior-information-enhanced Dynamic Imaging using Single or Multiple Coils 
with k-t BLAST and k-t SENSE. Proc Int Soc Magn Reson Med, Abstract 2369. 
2002). 

t 

5 

In the following, the k-t BLAST method is briefly reviewed, k-t SENSE is simiar, 
but it additionally takes the complex sensitivity information of the receiver colls into 
account, and it is described in detail in the aforementioned document. 

to Image content in a cine series of the heart exhibits a high degree of correiatior in 
both space and time. This is demonstrated in Figures 1 and 2. Figure 1 shows tie 
sampling pattern in k-t space {k = wavenumber, t = time), indicating that fr-space is 
fully sampled at every time point t The density and extent of sampling in /£space 
are chosen according to the desired field of view and resolution, respectively. 

.15 Figure 2 shows the corresponding fully-sampled, cine short-axis acquisition: of 1 ne 
heart 1. Image information along a vertical profile 2 through the heart (Figurie 2) is 
confined to compact regions in x-f space 3 (x - spatial position, f = temporal 
frequency), according to the dynamics of the corresponding image pixels. In : the k- 
t BLAST method, acquisition efficiency is increased by applying sparse sampli ig 

20 to k-t space according to a sheared grid pattern 5 (Figure 3c). This resujts in a 
tighter packing of signal replica in the reciprocal x-f space 6. Using low-resolution 
image information sampled for all cardiac phases 7 (Figure 3a), an estimate 8 of 
the expected signal distribution in x-f space is obtained, which is us ! ed to 
determine a reconstruction filter 9 (Figure 3b). The sparse sampling in the actual, 

25 high-resolution acquisition results in multiple replica 6 in x-f space, which fn turn 
manifest as multiple -fold-over in the actual image series 10 (Figure 3c) 
multiplying the acquired, high-resolution data 6 with the reconstruction filter: g, 
unfolded image series 1 1 is obtained after Fourier transform along the time 
(Figure 3d). 



The following describes the data acquisition scheme for a typical application jof the 
present invention: 3D single-breath cine SSFP cardiac imaging. In additioji, the 
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following describes the problem with a straightforward implementation 
describes the solution to overcome this problem. The solution is parti of 
present invention. . 



and 
the 



Samphng along a sheared grid pattern was implemented into the acquisition 
software of a Philips Gyroscan Intera system (Philips Medical Systems, Best, The 
Netherlands). An elliptical /c-space shutter 12 was applied to reduce the nurnbe of 
points to be sampled by approximately 25% by skipping the corners of 'tUpt v* 
(F«gure 4a). The straightforward approach of sampling *,space is to use i linlar 
profile order, which is the default order, but it leads to image artifacts 13 >igL 
4c). These artifacts originate from the fact that at the end of each ^ line 14 1 thj ^ 
position is incremented by one, but the k y position now jumps to the sWL 
position on the next k* line 15. This will be referred to as a line jump hereafter. 
Usually, when the. jump is small, the rewinding gradient is similar in magnitude but 
opposite in polarity to the phase-encoding gradient, so eddy currents generated by 
the former are well compensated by those of the latter. In the case of line Jumps, 
the eddy currents are not compensated, resulting in image artifacts. ThS sane 
also applies at the end of each cardiac phase 16; At that point, both the ky and k* 
positions jump back to the starting positions in the next cardiac phase if Tie 
jump is usually small for k* but can be quite large in k y , potentially leading to ed iy 
current problems. : ! 

For a well tuned gradient system, these artifacts should be relatively ^ms ill. 
Particularly for 2D imaging, these artifacts are almost at a negligible level g 
(Figure 5a) (see B.A Jung et. al. Magn. Reson. Med. 2002; 48, p. 921-92H). 
However, these artifacts are drastically amplified in case of accelerated acquisition 
20, particularly at high acceleration factors. This is due to two reasons. Firstly, k- 
space is sampled sparsely in order to accelerate the acquisition. This means thjat 
larger, jumps 21 are involved going from one profile to the next, thereby resul 
more residual eddy currents. Secondly, since fewer profiles are acquired 



ting 



in 



each 



time point or cardiac phase, the line jumps occur more frequently. For exdmpia, 
suppose that in real-time 2D imaging, N profiles are acquired. If the usuaJ linear 
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profile order is used, a Pine jump would occur after every N lines. If the 
is accelerated by, for example, 8-fold, a line jump would occur after every 
lines. Therefore, the large jumps in fr-space occur more frequently. .ajnd 
generated eddy currents remain significant for the most part of data acquisition 
the acquired data are contaminated strongly by image artifacts. : j 



acquisition 
V8 
the 
so 



To reduce eddy-current related distortions induced by the changing I phase 
encoding, gradients, the present invention employs a Known eddy ' current 
compensation technique (see B.A. Jung et. al. Magn. Reson. Med. 2002;! 48, p. 
921-925), Previously, this technique was only used for unaccelerated scans, tor 
which the eddy current artifacts remain relatively mild for a well tuned hakiwire 
system, as mentioned. Therefore, it was even mentioned that slich 
compensation technique was unnecessary in the case of 2D imaging (see B A. 



Jung et. al. Magn. Reson. Med. 2002; 48, p. 921-925). 



The present invention is further based on the observation that the eddy^cummt 
artifacts, are drastically amplified with accelerated data acquisition. As a result, an 
eddy-current reduction technique, such as the said technique, must be usejdj, evsn 
in the case of 2D imaging. In the said technique for eddy-current reduction, an 
alternating profile order is used to minimize the large jumps in /r-space. For 
example, in 2D real-time imaging, the profiles are ordered such that /r-space is 
swept along one direction at one time point, and along the opposite direction in 
the next time point (see arrow indicating seep direction In Figures 5c and 5oj. This 
profile order significantly reduces the image artifacts both in normal unaccelerated 
acquisition 22 (Figure 5c) and in highly accelerated (e.g. 8-fold) acquisftipn ;>3 
(Figure 5d). In cine 3D acquisition, consecutive profiles in the ky-ke plarle are 
traversed in a contiguous fashion with the overall direction of sampling being 
reversed for odd cardiac phases (Figure 4b). This alternating profile order sc heme 
for subsequent k^ profiles and cardiac phases reduces eddy-current [related 
distortions by ensuring maximum compensation of phase-encode gradient obes, 
which leads to substantial reduction in image artifacts (Figure 4d) 18. 



i 
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It should be noted that the large jumps in /(-space and the resulting eddyjcunent 
problems may also affect other Imaging applications as well, including! banille! 
imaging at high acceleration factors, and reduced field-of-view imaging 'with a 
highly reduced field of view. Previously, this problem had not been rcppcjrtec in 

5 those applications because they are generally unable to achieve the high 
acceleration factors used in the present invention. For parallel Imaging, this 
problem may begin to arise due to the recent interest in using higher magnetic 
field strengths. At higher field strengths, the achievable acceleration factor i 
increased, so the aforementioned problem is expected, and the proposed splut 

10 can be applied to mat application as well 

I 
i 

I 

Typical embodiment: 3D cine cardiac imaging : ! 

Using 5-fold undersampling (5x k-t BLAST), volumetric data sets covering tie 
heart with 20 slices at a spatial resolution of 2x2x5 mm 3 were recorded viith 20 
cardiac phases in a single breathhold of 20-22 sec. with the training data acquired 
in a separate 5 sec prescan. When the training data were acquired together with 
the undersampled data in a single scan (referred to as interleaved acquisition 
hereafter), the breathhold duration was prolonged slightly to 25-27 se<j. Tie 
central, densely-sampled region in the kyk* plane consisted of 49 frrofilas 
providing low-resolution images without aliasing, A schematic of the s&mpling 



pattern applied is given in Figure 6. The central portion 15 of the ky-k* sf^ce 



16 



;s 
is 



sampled at the full FOV in each cardiac phase whereas the outer portion 
under-sampled by a factor of five with the pattern being shifted as a function of 
25 cardiac phase (Figure 6b). The gray dots 24 indicate the sampled pointe, the 
dashed dots 25 represent points sampled in the previous cardiac iphasp. 
Remaining acquisition parameters were as follows: FOV: 320x210 mm 2 ,: jTE/TR: 
1.55/3.1 ms, flip angle: 45 deg M heart phase interval: 28-35 ms. Taking the 
acquisition of low-resolution training data into account, the net acceleration factor 

30 was 4.3. i j 

i 
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Navigator-based volume tracking (see e.g. M.V. McConnell et al. Prospecive 
adaptive navigator correction for breath-hold MR coronary angiography.! Magn. 
Reson. Med. 1997; 37, p. 148-152) was implemented in order to ensure 
consistency between the low-resolution, fully sampled data and the high- 
resolution, under-sampled data when the two datasets were acquired in separate 
acquisitions. Using pencil-beam excitation on the diaphragm; the navigator signal 
corresponding to the breath-hold position in the first acquisition was recorded smd 
correlated with a navigator acquired prior to the second acquisition. \^oluiT»e 
tracking was applied to account for possible shifts in breath-hold positions, j 

Imaoe reconstruction 



The data were divided into two sets for reconstruction: a training set (Figure 3a) 
and an under-sampled set (Figure 3c). The training set consisted of the dens sly 
15 sampled central region of fr-space only 4; while the under-sampled set consisted 
of all the data sampled along the sheared grid pattern 5. Both data seti were 
inverse Fourier transformed along all axes. The training set yielded a self of low- 



set 



of 



resolution images in x-f space 8, while the under-sampled set yielded a 
aliased images in x-f space 6. For each set of aliased voxels in x-f space, the 
20 aliasing was resolved by the following equation: 

p - p DC + |©E w (e©E H [a -a DC j 

Vector p represents the reconstructed (unallased) intensities. © is a diagonal 
matrix, with the diagonal elements representing the squared expected intensitie s, 

25 as determined from the training data. E denotes the encoding matrix, represent! lg 
the process of aliasing. ¥ is the noise variance of the acquisition data. Vejcto: 
holds the aliased intensity of the voxels. The time-averaged signals of the.jobject 
(a D c) were subtracted prior to reconstruction. After reconstruction, the; jobjc et 
corresponding to the time-invariant signals (ppc) was added back. The sijuare- 

30 bracketed portion corresponds to the reconstruction fitter 9 as shown in Figure 3I> 
The training data were substituted into the reconstructed image for improved data 
consistency, if the training data were acquired in an interleaved fashion. To reduce 

i 
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computation, reconstruction was performed separately for each coil element, und 
the reconstructed images were combined afterwards as the root-mean-s^uaro of 
the coil [mages. In general, k-t SENSE, which takes into account the cemllex 
sensitivities of the receiver coils, may also be used, instead of the coil-bylcoil k-t 
BLAST reconstruction described above. 

i 

|n-vivo validation . 

The study population consisted of four healthy volunteers and four jcarc iac 
patients undergoing cardiac MR imaging for assessment of congenial hoart 
disease. In both groups, conventional, multi-slice 2D SSFP acquisition's were 
obtained additionally for comparison, which required 9 sec per slice. Each s ice 
was acquired in a separate single breath-hold with Imaging parameters 
comparable to the single breath-hold 3D k-t BLAST acquisition. Left vehicular 
contours were extracted using standard post-processing software, and the fcontour 
area from the single-slice and 3D acquisitions was compared. 



Results 



3D data sets were successfully acquired in all subjects. Nearly araract-free 3D 
volumes were reconstructed, allowing reformatting along arbitrary views! Si ght 
image ghosting remained from residual eddy-current effects. In a few subjects with 
large variations of the R-R interval during the breath-hold, early diastolic (frames 
exhibited a small degree of blurring. 

i 

• i 

In Figure 7a, a cine 3D single breath-hold data set of central fifteen slice^ at ithe 
end systolic phase of. a 5x k-t BLAST acquisition obtained in a healthy subject is 
shown. Figures 7b and c depict selected time frames of a mid-ventricular sh art 
axis slice (Figure 7b) along with reformatted four-chamber views (Figure 75c). "he 
actual temporal resolution was 35 msec. 

' i 

■ I 

■ li 
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Rgure. 8 shows example Images obtained from one of the patients who had 
repaired Tetralogy of Fallot and pulmonary regurgitation causing right ventricular 
. dilation. In all patients, training data were obtained in a separate prescan, listing 5 
sec, in order to allow shorter breath-hold durations of the high-resolution 3D 
5 acquisition. The data set was obtained similar to the set of Fig. 7 with 5x k-t 

BLAST. The actual temporal resolution was 28 msec. i 

i 
i 

In in-vivo validation, the mean difference of left ventricular contours extracted in 
three exemplary slices of the 3D data set and the multi-slice 2D acquisition v/as 
10 3.41+6.52%. \ 

In the present invention it is shown that k-t BLAST can be applied successfully to 
cine 3D acquisition of the heart at 5-fold acceleration. The increased sc&n 
efficiency allows the acquisition of volumetric data sets of the heart at hlghi sps tial 
is (2x2x5 mm 3 ) and temporal (28-35 ms) resolutions in a single breath-hold! las ing 
20-22 sec (for a heart rate of 55-80 beats/min), with the training date? being 
obtained separately. The acquisition strategy is well tolerated by healthy sjubjejcts 
and patients. . 

20 The modified profile order according to the invention proved to be important euen 
for 2D imaging, particularly when high acceleration factors are used. The invention 
proved effective in reducing eddy-current related image artifacts by eliminating the 
frequent large jumps in fc-space associated with sparse sampling. In cine ufiaging, 
residual artifacts are considered to be related to the varying number of pummy 

25 cycles at the end of each cardiac cycle, depending on the R-R interval.' ^urtier 
improvements in eddy current reduction can be achieved by using dumrnyj eye les 
after each trigger to provide sufficient time for the eddy currents to decajr, or to 
smooth any fc-space jumps. ; 

i 
i 

30 To ensure consistency between breathhold levels when training data are ^bquired 
in a separate prescan, navigator-based volume tracking Is applied. The colreclion 
factor is set to 0.8, which is considered only a rough estimate of the cxyvelalion 
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between diaphragm and heart displacement. Inaccuracies, however, do not affect 
reconstruction to a large extend as training data are of very low spatial resolution 
(2x23x11 mm 3 ). ; 



Fig. g shows diagrammatically a magnetic resonance imaging system in wriich 



invention is used. 



the 



20. transmission and receiving coils 13, 15 for generating RF excitation pulses land 
picking up the magnetic resonance signals, respectively. The transmission coil 



The magnetic resonance imaging system includes a set of main coils 10 viheniby 
a steady, uniform magnetic field is generated. The main coils are constructed, for 
example in such a manner that they enclose a tunnel-shaped examination j space. 
The patient to be examined is slid on a table into this tunnel-shaped examination 
space. The magnetic resonance imaging system also includes a number kf gradi- 
ent coils 11,12 whereby magnetic fields exhibiting spatial variations, notably in the 
form of temporary gradients in individual directions, are generated so as to be 
superposed on the uniform magnetic field. The gradient coils 11, i!2 are 
connected to a controllable power supply unit 21. The gradient coils 11, !12 are 
energized by application of an electric current by means of the power supply unit 
21. The strength, direction and duration of the gradients are controlled by control 
of the power supply unit. The magnetic resonance imaging system alsoiirpcluiles 

for 
13 

is preferably constructed as a body coil whereby (a part of) the object to be 
examined can be enclosed. The body coil is usually arranged in the m^gn stic 
resonance imaging system in such a manner that the patient 30 to be examined, 
being arranged in the magnetic resonance imaging system, is enclosed j by the 
v body coil 13. The body coil 13 acts as a transmission aerial for the transmission of 
the RF excitation pulses and RF refocusing pulses. Preferably, the body coll 13 
involves a spatially uniform intensity distribution of the transmitted RF pulses, "he 
receiving coils 15 are preferably surface coils 15 which are arranged on! or nsar 
the body of the patient 30 to be examined. Such surface coils 15 have I a high 
sensitivity for the reception of magnetic resonance signals which is also isjpati ally 
inhomogeneous. This means that individual surface coils 15 are mainly sensrjve 
for magnetic resonance signals originating from separate directions, ij.e. fiom 
separate, parts in space of the body of the patient to be examined. The coil se|isi- 
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tivlty profile represents the spatial sensitivity of the set of surface coils. The trans- 
mission coils, notably surface coils, are connected to a demodulator 24 ind the 
received magnetic resonance signals (MS) are demodulated by means j of the 
demodulator 24. The demodulated magnetic resonance signals (DMS) are applied 
5 to a reconstruction unit. The reconstruction unit reconstructs the magnetic reso- 

* i 

nance image from the demodulated magnetic resonance signals (DMp) sind 
optionally on the basis of the coil sensitivity profile of the set of surface coils. 1 he 
coil sensitivity profile has been measured in advance and is stored, for example 
electronically, in a memory unit which is included in the reconstruction Unit. The 

io reconstruction unit derives one or more image signals from the demodulated 
magnetic resonance signals (DMS), which image signals represent one of mc re, 
possibly successive magnetic resonance images. This means that thej signal 
levels of the image signal of such a magnetic resonance image represent ;he 
brightness values of the relevant magnetic resonance Image. The reconstruction 

is unit 25 in practice is preferably constructed as a digital image processing unit 25 
which is programmed so as to reconstruct the magnetic resonance image fijom he 
demodulated magnetic resonance signals and optionally on the basis of the ;oil 
sensitivity profile. The digital image processing unit 25 is notably programmed 



so 



as to execute the reconstruction in conformity with the present invention. " "he 
20 image signal from the reconstruction unit is applied to a monitor 26 so that the 
monitor can display the image information of the magnetic resonance! image 
(images). It is also possible to store the image signal in a buffer unit 27 wpile 

awaiting further processing, for example printing in the form of a hard copy.: 

i 
i 

25 In order to form a magnetic resonance image or a series of successive magnetic 
resonance images of the patient to be examined, the body of the pajtienv is 
exposed to the magnetic field prevailing in the examination space. The ^teady, 
uniform magnetic field P i.e. the main field, orients a small excess number of the 
spins in the body of the patient to be examined in the direction of the main fiiild. 

30 This generates a (small) net macroscopic magnetization in the body. Thesb sp ins 
are, for example nuclear spins such as of the hydrogen nuclei (protons), but e 
tron spins may also be concerned. The magnetization is locally influenced 
Explication of the gradient fields. For example, the gradient coils 12 apply ; a se 

tion gradient in order to select a more or less thin slice of the body. Subsequently, 

i i 
• i 

. i 

! I 
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the transmission coils apply the RF excitation pulse to the examination s^ace in 
which the part to be imaged of the patient to be examined is situated, hjhe RF 
excitation pulse excites the spins in the selected slice, i.e. the net magn^feat on 
then performs a processional motion about the direction of the main field, bur ng 
this operation those spins are excited which have a Larmor frequency within 



:he 
so 



frequency band of the RF excitation pulse in the main field. However, it jis a 
very well possible to excite the spins in a part of the body which is much laroer 
man such a thin slice; for example, the spins can be excited in a thre4-dim *n- 
sional part which extends substantially in three directions in the body. After the RF 
excitation, the spins slowly return to their initial state and the macroscopic m ag 
netization returns to its (thermal) state of equilibrium. The relaxing spins tften emit 
magnetic resonance signals. Because of the application of a read-out gradient 
and a phase encoding gradient, the magnetic resonance signals have a plurality of 
frequency components which encode the spatial positions in, for example he 
selected slice. The fr-space is scanned by the magnetic resonance signals by 
application of the read-out gradients and the phase encoding gradients. According 
to the invention, the application of notably the phase encoding gradients rekuitn in 
the sub-sampling of the fc-space, relative to a predetermined spatial resoiiitior of 
the magnetic resonance image. For example, a number of lines which is fp\> small 
20 for the predetermined resolution of the magnetic resonance image, for example 
only half the number of lines, is scanned in the fc-space. ! 
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Claims 

1. A magnetic resonance imaging method to produce successive 
resonance images wherein: 



magnetic 



a series of. successive magnetic resonance signals is obtained by steajdy- 
state free precession imaging, 

successive sets of magnetic resonance signals are acquired by successively 
scanning respective sets of points in K-space in an undersampled fasliion, 
the magnetic resonance signals are acquired in conjunction with art edjdy- 
current reduction technique, and j 
successive magnetic resonance images are reconstructed frbjm 



successive sets of magnetic resonance 
reconstruction method. 



using a 



the 
Suitable 



A magnetic resonance imaging method as claimed In Claim 1, wherein 
eddy-current reduction technique employs alternating sweep directions 
sampling k-space. ; 



he 
in 



A magnetic resonance imaging method as claimed in Claims 1 or 2, where in 
successive sets of magnetic resonance signals are acquired by successively 
scanning respective sets of points in k-space in an undersampled: fashion 
such that the ensemble of successive sets cover the entire .portion oi k- 
space at full sampling density , j j 

successive updates of a training set of magnetic resonance signals are 
obtained from the magnetic resonance signals, either in the same span o r in 
a separate scan, by further acquisition of the central portion of k-space at tull 
sampling density or with slight undersampling if multiple receiver antennae 
are. used, 

the undersampled sets of magnetic resonance signals are successively 
updated by further undersampled scans of the entire k-space, | j 



i 
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• i 

! i 



resonance 



data, 



a baseline Image Is optionally reconstructed from the training data 
under-sampled data, or from data acquired separately during timai jj 
with little or no motion, 

a distribution of likelihood of changes in the successive magnetic 
images is identified from the static reference image and/or the training 
in the space spanned by geometrical space alone or by geometrical 
and temporal frequency, • j 

successive magnetic resonance images are reconstructed frbjn 
respective sets of undersampled magnetic resonance signals on the .paste 
the identified distribution of likelihood of changes, and optionally the Basel 
image, and 

the magnetic resonance signals are optionally acquired by way of a 
antennae system having a spatial sensitivity profile, and the su 
magnetic resonance images are reconstructed from the respective: Lets 
undersampled magnetic resonance signals on the additional and ojp 
basis of the sensitivity profile of the receiver antennae. 



space 



the 
of 
ne 



receiver 
;eessiive 



of 



tional 



4. A magnetic resonance imaging method as claimed in Claims 1 or 2, whergin 
the magnetic resonance signals are acquired by way of a receiver ahj:e 
system having a spatial sensitivity profile, and successive mag. 
resonance images are reconstructed from the respective sets 
undersampled magnetic resonance signals on the basis of the sensitivity 
profile of the receiver antennae. 



t 

4 I 

A magnetic resonance imaging method as claimed in Claims 1 or 2, Uh 
successive magnetic resonance images are reconstructed from 
respective sets of undersampled magnetic resonance signals on the b'asis 
a reduced field of view, where changes in image contents are assumed 
take place. . 



A magnetic resonance imaging method as claimed in one of Claims |l 
wherein an elliptical k-space shutter is applied. j ; 



tc 
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7. A magnetic resonance imaging method as claimed in one of Claims 1 to 
wherein navigator-based volume tracking is applied. . 



6, 



8. A magnetic resonance imaging system arranged to 

obtain a series of subsequent magnetic resonance signals by steady-state 
free precession imaging; 

apply an eddy-current reduction technique, such as by alternating the swebp 
directions of sampling in k-space, ! 
acquire a set of magnetic resonance signals in an undersampled fashion, 
optionally acquire the magnetic resonance signals by way of a receiver 
antennae system having a spatial sensitivity profile, j 
optionally acquire an additional training set of magnetic resonance sigrials, 

1 • i 

optionally reconstruct a baseline image from the training data jand/or 
undersampled data, or from data aquired separately during time periods \A|ith 
little or no motion, j 
optionally identity a distribution of likelihood of changes in the success 
magnetic resonance images from the baseline image arid/or the .t^ain 
data, in the space spanned by geometrical space alone or by geonptrif^l 
space and temporal frequency, 

optionally reconstruct the successive magnetic resonance images fnbrn 1 he 
respective sets of magnetic resonance signals of the dynamic series jon the 
basis of: 

- the identified distribution of likelihood of changes J 

- the baseline image. j 

- the sensitivity profile of the receiver antennae, and/or . j 

- a reduced field of view where changes in image contents are assuhnedj to 

• ! 

take place. 



9. A computer programme product comprising instructions to « j 

obtain a series of subsequent magnetic resonance signals by steady-stiate 



ve 
ng 



free precession imaging, 
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apply an eddy-current reduction technique, such as by alternating the* swepp 
directions of sampling in k-space, i 
acquire a set of magnetic resonance signals in an undersampled fashion, 
optionally acquire the magnetic resonance signals by way of a fecei^er 
antennae system having a spatial sensitivity profile, : j 

optionally acquire an additional training set of magnetic resonance signals, 
optionally reconstruct a baseline image from the training data Widbr 
undersampled data, or from data aquired separately during time periods wjith 
little or np motion, 1 j 

optionally identity a distribution of likelihood of changes in the successive 
magnetic resonance images from the baseline image and/or the training 
data, in the space spanned by geometrical space alone or by geometrical 
space and temporal frequency, j 
optionally reconstruct the successive magnetic resonance images fhim the 
respective sets of magnetic resonance signals of the dynamic series son the 
basis of: 

i 

- the identified distribution of likelihood of changes j 

t 

- the baseline image, 

- the sensitivity profile of the receiver antennae, and/or j 

i 

- a reduced field of view where changes in image contents are assujneq to 
take place. j 

i 
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Summary 

» 

This invention describes the combination of SSFP, a method for accelerating ds ta 
acquisition, and an eddy current compensation method. This synergistic 
combination allows acquisition of images with high signal-to-noise ratio, hijh 
image contrast, high spatial and temporal resolutions, and good immunity Against 
system instabilities, k-t BLAST and k-t SENSE are the preferred method or 
accelerating data acquisition, since they allow high acceleration factors, b<4 ott er 
methods such as parallel imaging and reduced field-of-vlew imaging arfe also 
applicable. Typical applications of this invention include cine 3D cardiac imjagir g, 
and 2D real-time cardiac imaging. j 



(Fig. 



3) 



• i 



i ! 



IBTp002 

Fmrvf TO < +;n9/f]7/90n3 19:15 



FmPf.nr.:578 P. 027 



02/07/2003 IB: 01 +41-56-4859495 



SPIERENBURG PARTNER 



1/6 



s + <p+B*on*»c*acc * ♦ » 6 » <> «> « 
♦ 9<>«tO a * ♦ j? * v ^ • * • « * 



O9»ac*-f*t"»o»<j * * * « » ♦ * n 
«A»«»4^«»4«Pf 9f m + mm • 
qi*»<f.o<>**»«iai * a * D • • o » » 



Fig. 1 




Fig. 2 



Fmpf.7eit:02/07/2003 19:16 



Empf.nr.:578 P. 028 



02/07/2003 18:01 +41-56-4853495 



SPIERENBURG PARTNER 




02/07/2003 18:01 +41-5G-4959495 



SPIERENBURQ PARTNER 




Fig. 6 



ii 

; \ 
' i 

!f 

• ■ 



02/87/2083 18:01 +41-56-4859495 



SPIERENBURG PARTNER 




This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



